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Papers  in  thla  series  have  resulted  from  technical  activities  of  The 
Hydrologic  Engineering  Center.  Versions  of  some  of  these  have  been  published 
in.  technical  journals  or  in  conference  proceedings.  The  purpose  of  this 
series  is  to  make  the  information  available  for  use  in  the  Center's  training 
program  and  for  distribution  within  the  Corps  of  Engineers. 
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ABSTRACT 

A  hydrologic-economic  simulation  nodal  waa  developed  to  avaluata  alter- 
nativa  protactlon  achaaaa  In  tha  daalgn  of  an  authorized  federal  flood 
control  project  for  125  alien  of  the  Tlbbee  River  flood  plain  in  Mlaelsalppl. 
The  aodel  requires  input  consisting  of  unit  hydrographs,  streanf low  routing 
coefficients  and  storage  functions,  a  pattern  atom,  rainfall  loss  rate 
functions,  and  flow-daaage-frequency  relations.  A  single  synthetic  pattern 
atom  was  used  in  conjunction  with  flow-frequency  curves  at  index  locations 
to  generate  a  series  of  floods  for  comparing  alternative  protection  schemes 
with  existing  conditions.  The  effect  of  channel  improvements  on  flood 
runoff  characteristics  was  evaluated  by  using  storage  routing  functions 
that  account  for  changes  in  storage-discharge  relations.  Based  on  results 
obtained  from  using  the  model,  a  channel  improvement  plan  was  tentatively 
selected  for  the  Tibbee  River  basin  from  the  alternative  schemes  evaluated. 
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INTRODUCTION 

A  federal  project  involving  extensive  stream  clearing, 
straightening  and  enlargement  was  authorized  by  Congress  in  1938. 

The  purpose  of  the  project  is  to  achieve  flood  control  and  drainage 
on  466  miles  of  flood  plain  along  22  tributary  streams  of  the  upper 
Tombigbee  River  in  Mississippi  and  Alabama.  A  current  study  conducted 
on  a  cooperative  basis  by  the  Mobile  District  and  The  Hydrologic 
Engineering  Center,  Corps  of  Engineers,  is  concerned  with  a  portion 
of  the  overall  project,  namely  125  miles  of  flood  control  improvements 
authorized  for  a  system  of  10  streams  in  the  Tlbbee  River  basin, 
Mississippi  (figure  l).  The  objectives  of  the  study  are  to  determine 
the  effect  of  alternative  channel  improvement  schemes  on  flood  runoff 
characteristics  in  the  stream  system,  to  determine  the  reductions  in 
flood  damage  associated  with  the  alternative  plans,  and  to  select 
from  among  the  alternatives  the  plan  that  would  most  economically 
achieve  the  flood  control  objectives. 

The  methodology  adopted  for  the  study  required  the  development 
of  a  generalized  hydrologic-economic  mathematical  model  to  simulate  the 
flood  runoff  characteristics  of  the  basin  under  existing  conditions  and 
for  alternative  improvement  schemas,  including  detention  structures 
that  are  planned  by  the  U.  S.  Soil  Conservation  Service .  This  paper 
describes  the  simulation  model  and  study  procedures  and  presents  soem 
of  the  results  of  the  study. 
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plain  of  the  10  streams  authorised  for  improvement .  The  basin  has 
a  maximum  width  of  30  miles,  an  average  width  of  21  miles,  and  a 
stream  length  of  7k  miles  (from  the  mouth  of  Tibbee  River  to  the  head 
of  Sakatonchee  Creek,  its  longest  tributary).  The  terrain  is  fairly 
rugged  with  rather  high  relief  varying  from  elevation  160  to  400  feet, 
m.s.l. 

The  Tibbee  River  is  formed  at  the  junction  of  Line  and  Sakatonchee 
Creeks  and  it  flows  easterly  about  2k  alias  to  join  the  Tombigbee  River 
through  the  right  bank  about  35  miles  northwest  of  Columbus,  Mississippi 
It  has  an  average  gradient  of  1.5  f««t  per  mile.  Its  tributaries 
range  in  length  from  about  5  miles  to  50  miles  with  gradients  generally 
steeper  than  that  of  the  main  stem.  The  flood  plain  along  the  Tibbee 
River  and  its  tributaries  averages  about  1  mile  in  width  aad  is  subject 
to  frequent  flooding. 

The  basin  area  is  located  within  the  Tomb i ghee  Hills  and  Black 
Prairie  districts  of  the  Gulf  Coastal  Plain  physiographic  province. 

The  downstreem  U  or  5  miles  of  the  Tibbee  River  traverses  the  Tombigbee 
Hills  district  with  topography  ranging  from  low,  smoothly  rounded 
hi  Lie  to  hi lie  and  ridge*  separated  by  narrow  valleys.  The  upstream 
reach  of  the  Tibbee  River  end  its  tributaries  is  located  in  the  Blaek 
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Prairie  district  where  the  topography  varies  fraa  nearly  flat 
to  — oothly  rounded  hi lie  or  low  relief. 

CUnatology.  The  Tibbee  River  baeln  has  a  temperature  eiiaate 
with  mmi  turners  and  ai Id  winters.  The  normal  annual  tenperature 
it  61s  degrees  with  nonthly  nonsals  ranging  free  kj  degrees  In 
January  to  Bl  degrees  in  July.  The  ainlsar  recorded  tenperature 
it  -10  degrees  and  the  naxlevss  is  113  degrees.  The  frost -free 
period  normally  lasts  from  April  to  Wovenber. 

Precipitation  during  the  year  is  abundant  and  fairly  well 
distributed.  The  normal  annual  rainfall  is  almost  49  inches,  of 
which  5B  percent  fails  in  the  winter  and  spring,  24  percent  in  the 
simmer,  and  18  percent  in  the  fall.  The  average  annual  snowfall  la 
about  35  inches. 

Flood  producing  storms  nay  occur  at  any  time  but  are  acre 
frequent  during  the  winter  and  spring.  Such  storms  are  usually  of  the 
frontal  type  covering  large  areas  and  lasting  from  2  to  4  days.  Summer 
storms  are  generally  of  the  thunderstorm  type  with  high  intensity 
over  email  areas.  Recorded  rainfalls  lit  the  basin  Include  lmximmi  of 
3.20  inches  in  a  1-hour  period,  6.53  inches  in  12  hours,  and  7.1»3 
inches  in  24  hours,  fle  l-  and  12-hour  maxisnass  occurred  during  a 
storm  on  June  1-2,  1947,  and  the  24 -hour  maximal  occurred  in 
January  1950. 

4*hd  Use .  About  half  of  the  total  basin  area  is  in  farms.  Sixty 
percent  of  the  farmland  acreage  is  about  equally  divided  between 
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cropland  and  other  uaaa  such  as  houat  lota,  roads,  pasturaa  and 
vaatelaad.  Tha  remaining  Uo  per cant  la  woodland.  Mora  than 
80,000  acre a  are  In  the  flood  plain  along  the  10  project  etreana . 

It  la  aatlnated  that  30  percent  of  thia  flood  plain  land  la  in 
cultivation,  1$  percent  ia  in  paature  or  idle,  and  55  percent  ia  wood¬ 
land.  It  ia  eatinated  that  flood  damages  to  flood  plain  development; 
average  $1,837,000  a  year,  consisting  mainly  of  agricultural  losses. 

The  improvement  of  the  project  streams  would  reduce  these  flood 
dasuges  and  permit  more  intensive  agricultural  use  of  the  flood 
plains. 

INVESTIGATION 

General.  To  accomplish  the  study  objectives,  it  was  concluded 
that  a  systems  analysis  of  various  alternative  schemes  of  basin  develop¬ 
ment  for  the  Tibbee  River  was  required,  and  that  a  hydrologic-economic 
simulation  model  should  be  developed  to  facilitate  the  analysis.  It  was 
recognised  that  extensive  channel  improvement  works  in  a  river  system 
could  cause  an  adverse  peaking  effect  on  runoff;  i.e.,  peak  rinoff  rates 
could  likely  be  larger  because  of  the  more  rapid  concentration  of  runoff 
resulting  from  more  efficient  conveyance.  Therefore,  the  model  was 
developed  primarily  for  the  purpose  of  evaluating  the  effects  of 
alternative  channel  Improvement  plans  on  the  flood  runoff  characteristics 
of  timing  and  magnitude  of  peak  runoff  rate.  Another  major  purpose  of  the 
model  was  to  determine  the  reductions  in  flood  damage  for  the  alternative 


plana.  A  generalised  computer  program  ^  (rainfall- runoff  nodal) 
developed  by  The  Hydrologic  Engineering  Cantor  wi  uaad  for 
modeling  purpoaaa. 

Tha  application  of  tha  gonoralisad  aiaulation  nodal  required 
tha  davalopnant  of  certain  hydrological  and  acononic  inputs,  such  as 
(l)  unit  hydrograph  and  loss  rata  criteria,  (2)  channel  routing  criteria, 

(3)  synthetic  pattern  atom,  (U)  di  jeharga- frequency  and  stage-discharge 
relations,  and  (5)  stage-daaage  relations. 

Economic  Studies.  Tha  procedure  uaad  for  data raining  reductions 
in  flood  daaage  attributable  to  change 1  improvement  works  was  to 
compute  tha  di f farence  in  awaraga  annual  flood  damages1  between 
existing  conditions  and  improved  conditions.  Average  annual  flood 
damages  for  existing  conditions  were  computed  from  known  or  derived 
stage-discharge,  discharge- frequency,  and  stage-damage  relationships. 

For  mproved  conditions,  average  annual  flood  damages  are  usually  computed 
from  these  saam  relationships  with  the  stage-discharge  relatlonehlp  modified 
to  reflect  the  changes  in  stage.  However,  for  this  study  it  was  also 
necessary  to  modify  the  discharge-frequency  relationship  to  reflect  the 
magnitude  of  the  peaking  effect  on  runoff  resulting  from  upstream  channel 
improvement  works.  The  need  for  modifying  this  relationship  made  the 
determination  of  flood  damage  reduction  more  difficult. 

In  formulating  the  optimal  plan  of  channel  iaproveamnt  works  for 
the  Tibbee  River  basin,  estimates  of  flood  dsmsgs  reduction  attributable 

1  Average  annual  flood  damage  ie  defined  as  tbs  expected  value  of  anmmsl  flood 
demage  computed  by  integrating  the  damage -ex oasis ace  frequency  function. 


to  each  alternative  plan  evaluated  were  determined.  For  this 
determination,  the  flood  plains  of  the  Tibbee  River  and  its 
tributaries  were  subdivided  into  21  planning  reaches  varying  in  length 
from  2  to  17  miles  (figure  l).  Based  on  detailed  economic  studies, 
flood  damages  were  determined  for  each  planning  reach  and  related 
to  stage  at  the  index  station  (representative  reference  point) 
selected  for  the  reach.  Using  stage-discharge  and  discharge-frequency 
relationships,  average  annual  flood  damages  were  computed  at  each 
index  station  for  existing  conditions  and  for  each  alternative  plan 
as  described  later  in  the  paper. 

Hydrologic  and  Hydraulic  Studies.  A  discharge-frequency 
statistical  analysis  was  made  from  streamflov  data  for  gaging  stations 
located  on  the  Tibbee  River  and  its  tributaries  and  on  adjacent  streams 
of  similar  size  sod  comparable  hydrologic  characteristics.  A  relationship 
of  drainage  area  versus  the  annual  mean  was  developed.  A  skew  coefficient 
of  zero  and  a  standard  deviation  of  0.2b  were  adopted  and  used  with  this 
relationship  for  determining  the  discharge-frequency  relationship  at 
each  index  station  for  existing  conditions  in  the  basin. 

Extensive  backwater  computations  for  existing  conditions  and  for 
alternative  plans  of  improvement  were  made  to  determine  stage-discharge 
relations  at  each  index  station,  pertinent  water  surfaoe  profiles,  and 
storage-outflow  relations  for  the  selected  channel  routing  reaches.  Data 
on  surveyed  cross -sections  aad  topographic  aape  were  available.  The 
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roughness  coefficients  used  in  Manning's  formula  were  based  on 
known  high  water  marks  and  rating  stations  in  the  basin. 

For  study  purposes ,  the  Tibbee  River  basin  was  subdivided  into 
216  subareas.  Synthetic  unit  hydrographs  were  derived  for  each 
subarea  by  the  Clark  unit  hydrograph  method  ^  using  generalised 
coefficients.  Unit  hydrograph  and  rainfall  loss  rate  coefficients 
were  derived  for  22  stream  gaging  locations  in  the  Tibbee  River  basin 
and  other  adjacent  stream  basins  from  an  analysis  of  streamflow  and 
rainfall  data  for  a  total  of  89  runoff  events  (2  to  8  per  location). 

The  results  of  this  analysis  were  used  for  developing  generalized 
coefficients  for  application  to  ungaged  locations. 

The  methodology  adopted  for  evaluating  the  effect  of  channel 
improvements  on  flood  runoff  characteristics  required  the  development 
of  a  single  representative  synthetic  pattern  storm  for  the  Tibbee  River 
basin.  It  also  provided  for  computing  a  synthetic  pattern  flood  for 
each  subarea  from  the  pattern  storm  and,  using  multiples  of  the 
pattern  flood,  for  determining  other  floods  of  pertinent  magnitudes  or 
frequencies. 

A  uniform  areal  distribution  of  rainfall  was  used  for  computing 
the  pattern  storm.  In  this  way,  each  portion  of  the  basin  had  a 
representative  effect  on  flood  determinations.  Although  a  single  pattern 
storm  nay  not  be  totally  representative  for  the  entire  Tibbee  diver 
basin  as  well  as  the  headwater  areas,  it  would  haw  been  extremely 
difficult  end  laborious  to  attempt  many  representative  storm  centerings 
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within  the  basin.  Also,  it  is  doubtful  that  results  would  have 
been  substantially  different  than  those  obtained  from  a  single 
uniform  storm  pattern.  While  the  effects  of  tributary  runoff  would 
actually  va~y,  depending  on  storm  centering,  the  net  effect  of  all 
centerings  and  all  storm  magnitudes  would  probably  not  be  substantially 
different  from  that  computed  using  a  uniform  distribution. 

The  area-depth-duration  relation  adopted  for  the  pattern  storm 
was  of  particular  importance  because  the  resulting  pattern  floods 
were  to  be  representative  of  historical  floods.  It  was  determined 
that  if  the  peak-to  volume  relationships  of  the  pattern  floods  are 
representative  of  historical  floods  throughout  the  basin,  then  other 
floods  derived  as  multiples  of  the  pattern  floods  would  be  equally 
representative  because  peak-to-volume  ratios  of  historical  floods  are 
relatively  independent  of  flood  magnitude.  The  adopted  criteria  for  the 
storm  were  a  rainfall  intensity  of  5-year  recurrence  interval,  a  rain¬ 
fall  duration  of  48  hours,  and  an  area  size  of  200  square  miles.  Based 
on  rainfall  data  in  Weather  Bureau  Technical  Papers  4o  and  U9,  it  was 
determined  that  the  total  storm  rainfall  was  5-7  inches. 

For  this  study,  the  pattern  flood  at  any  location  in  the  basin  is 
defined  as  the  runoff  from  the  pattern  storm.  The  pattern  flood 
hydrograph  for  each  subarea  was  computed  by  applying  the  pattern  storm 
rainfall  excess  to  the  synthetic  unit  hydrograph.  The  pattern  flood 
hydrographs  for  the  other  locations  in  the  basin  were  determined  by 


routing  an  i  r<jmturiing  these  suterea  pattern  flood  hydrographs.  The 
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Muskingum  routing  method  was  used  for  the  smaller  streams  in  the  head¬ 
water  areas  and  the  nonlinear  storage  routing  method  (modified  Puls) 
was  used  for  the  other  stream  reaches.  The  latter  method  required  the 
development  of  storage-outflow  relations  for  each  selected  channel 
routing  reach  for  existing  and  improved  conditions.  In  this  relation, 
the  storage  value  is  the  volume  of  water  in  the  channel  and  overbank 
within  the  routing  reach  corresponding  to  a  given  steady  state  water 
surface  profile.  The  storage-outflow  relation  for  existing  conditions 
was  adjusted  to  reflect  the  changes  in  storage  for  each  plan  of 
channel  improvement  considered.  Therefore,  although  the  same  pattern 
storm  was  used,  the  resulting  pattern  floodB  at  a  given  location  were 
different  for  each  alternative  plan  as  illustrated  in  figure  2  for 
one  of  the  index  stations  on  Catalpa  Creek. 

A  series  or  set  of  nine  synthetic  flood  hydrographs  at  all 
locations  in  the  basin  was  developed  for  existing  conditions  and  for 
each  considered  plan  of  improvement.  The  set  of  flood  hydrographs  was 
synthesized  by  multiplying  all  of  the  subarea  pattern  flood  hydrograph 
ordinates  by  nine  predetermined  ratios  and  by  routing  and  combining 
the  subarea  pattern  floods  computed  for  each  multiple.  The  same 
multiples  used  for  existing  conditions  were  used  for  each  plan  of 
Improvement.  Multiples  were  selected  to  produce  flood  magnitudes 
corresponding  to  a  frequency  range  from  four  to  twelve  times  per  year 
(zero  flood  damage  stage)  to  approximately  once  in  200  years. 
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CONCLUSIONS 


The  Boat  important  concept  of  tha  investigation  waa  the 
determination  of  tha  runoff  peaking  affaet  of  channal  iaprovement 
work:*  baaed  on  changes  in  tha  storage -outflow  ralation  of  tha  flood 
plain  and  atraaaa.  Anothar  important  and  uniqua  concapt  waa  tha 
computation  of  aynthatic  flood  hydrographs  at  nuaaroua  locationa  aa 
Bultlplaa  of  pattam  flood*  derived  from  ona  single  rapraaantatlve 
aynthatic  pattam  atom  for  a  atraan  ay  a  tan. 

Uoatraam  channal  inprovanant  worka  for  flood  control  would  tand  to 
incraaaa  tha  magnitude  of  paak  flow  rata*  now  axparlancad  on  tha  Tibbea 
Hi var  main  stem.  Thay  would  alao  tand  to  incraaaa  tha  average  annual 
flood  dasmges  on  tha  main  atan  unlaaa  it  waa  enlarged  to  accoaamdata  the 
increased  discharges.  Tha  aagnitude  of  thaaa  ef facta  waa  aaaaaaad  with  a 
raaaonabla  Mount  of  computation  using  a  aiaplifiad  hydro logic -a eoooaic 
ainulation  modal  nf  tha  baain,  along  with  certain  neceaaary  aaauaptiona 
and  generalisations. 

The  procedure  deacribed  herein  for  tha  evaluation  -f  flood  control 
channels  in  rural  areas  could  alao  ha  appliad  to  tha  study  of  urban  flood 
problems.  The  genera l izat ions  would,  of  course,  ha  tailored  to  the 
specific  area  of  study,  but  tha  methodology  and  computational  procedure 
would  be  the  s«uae.  Future  studies  should  investigate  the  possibility  of 
adding  to  tha  modal  an  optimisation  routine  to  provide  for  the  selection 
of  projects  involving  a  cos&lnation  of  two  or  aora  structural  measures 
(dentention  structures  together  with  levees  or  channels,  ate.)  which  are 
to  achieve  a  given  objective  (protection  against  a  specif iad  flood). 
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